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HISTORICAL  DEVELOPMENT 
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The  REEDM  program  has  its  roots  in  modeling  efforts  that  go  back  more 
than  20  years.  Beginning  with  a  beryllium  hazard  study  for  the  Pacific 
Missile  Range  (PMR)  in  1963  and  continuing  with  the  current  effort  to 
predict  acid  drop  deposition  from  Shuttle  operations  in  the  relatively 
complex  terrain  at  VAFB.  As  early  as  the  late  1960's,  it  was  recognized 
that  existing  point  and  volume  source  diffusion  models  were  inadequate  to 
handle  the  large  spatial  scales  encountered  in  the  launch  of  space 
vehicles.  By  1970  the  multilayer  model  had  been  formulated  and  implemented 
in  a  computer  program  to  be  used  by  NASA.  One  of  the  first  uses  of  the 
model  was  the  assessment  of  the  Titan  III  D  toxicity  hazard  at  VAFB.  The 
first  version  of  the  model  was  very  primitive,  requiring  specification  of 
stabilized  cloud  parameters  and  direct  input  of  meteorological  data.  With 
time  a  preprocessor  was  developed  which  calculated  the  cloud  height  as  a 
function  of  time,  the  stabilized  source  dimensions,  turbulence  parameters, 
and  wind  profiles.  The  model  now  Included  the  capability  to  handle 
abnormal  launches  and  had  stored  parameters  to  describe  Titan,  Delta-Thor, 
Mlnuteman,  and  Space  Shuttle  launches. 


By  the  end  of  1982  the  REEDM  program  was  operational  at  KSC  and  had 
supported  launches  of  both  the  Titan  and  the  Space  Shuttle.  A  user's 
manual  had  been  written  and  most  program  inputs  were  either  selected 
according  to  the  launch  being  supported  or  were  available  from  Information 
on  file.  A  program  execution  could  be  run  for  HCl  concentration 
calculations  with  every  parameter  defaulted  except  for  the  rawlnsonde  input 
file  and  the  height  of  the  mixing  layer  and  be  expected  to  produce  correct 
answers. 


At  this  time  it  was  recognized  that  the  fallout  of  acid  drops  from  the 
ground  cloud  was  a  major  problem;  one  that  could  not  be  addressiid  by  the 
current  version  of  REEDM.  In  1983  and  '84  algorithms  to  predict  the 
deposition  of  the  acid  drops  were  developed  and  Incorporated  into  the 

HP-1000  mini-computer  version  of  REEDM.  «  -  -  _ 
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2. 


PROGRAM  COMPONENTS 


The  REEDM  program  is  logically  divided  into  5  parts:  meteorological 
inputs,  source  inputs  dependent  on  launch  vehicle  and  type  of  launch, 
cloud-rise  and  material  distribution  algorithms,  the  dispersion  model 
algorithms  (there  are  three  -  -  dosage/concentration,  gravitational 
deposition,  and  washout  deposition),  and  output  routines. 

2.1  Meteorological  inputs 

The  program  meteorological  inputs  all  come  from  disk  files  resident  on 
the  same  computer  system  as  the  program.  During  operational  launch  support 
the  meteorological  data  is  updated  by  other  support  programs  which  collect 
the  data  from  rawinsonde  ascents,  tower  data,  and  accoustic  radar.  Not  all 
data  sources  are  currently  available  at  KSC.  Although  it  is  desireable  to 
use  all  the  sources  of  meteorological  data  listed  above,  the  REEDM  program 
has  been  developed  to  execute  with  only  the  rawinsonde  in  order  to  continue 
to  provide  output  for  hazard  predictions.  The  capability  at  VAFB  is 
increased  by  the  addition  of  a  mesoscale  wlndfleld  model  which  will 
incorporate  the  Influence  of  the  complex  terrain  found  there. 

2.2  Source  Inputs 

Source  Inputs  are  selected  from  stored  values  by  the  choice  of  launch 
vehicle  and  launch  conditions  and  consist  of  trajectory  data,  heat  emission 
data,  and  exhaust  chemical  constituents.  These  data  affect  the  cloud-rise 
and  total  quantity  of  each  reaction  product  used  as  input  to  the  dispersion 
models. 

2.3  Cloud  rise  and  other  vertical  variable  algorithms 

The  cloud  rise  algorithm  is  based  on  the  work  of  Briggs  and  for  normal 
launches  uses  his  Instantaneous  cloud  rise  model.  This  model  assumes  that 
entrained  air  increases  the  radius  of  the  ascending  cloud  as  a  linear 
function  of  the  height  gained  since  cloud  formation.  The  REEDM  program 


Is  a  multi-layer  model  with  the  layer  boundaries  defined  by  the  rawinsonde 
reporting  levels  entered  in  the  data  file.  In  order  to  reduce  the  error 
which  would  occur  with  excessively  thick  layers,  the  program  will 
interpolate  intermediate  levels  as  needed.  During  incorporation  of  the 
acid  drop  deposition  algorithm  into  REEDM,  we  found  that  limitations 
also  had  to  be  placed  on  the  amount  of  directional  shear  within  a  layer. 
Once  a  satisfactory  layering  has  been  found  by  the  program,  the  cloud  rise 
algorithm  computes  the  cloud  height  and  vertical  velocity  as  a  function  of 
time  with  output  values  at  each  layer  boundary  below  the  cloud 
stabilization  height.  The  calculation  of  cloud  vertical  velocity  is  done 
on  a  layer-by-  layer  basis  with  the  value  of  the  stability  parameter  used 
in  the  cloud  rise  algorithm  for  each  layer  being  computed  from  the  average 
of  the  stability  parameter  through  the  vertical  extent  of  the  cloud  (from 
the  bottom  of  the  cloud  through  the  top  of  the  cloud).  The  heat  available 
for  cloud  rise  through  a  layer  is  computed  from  the  stored  power  law 
approximation  to  the  launch  vehicle  ascent  height  vs  time  and  the  heat 
emission  rate.  The  heat  available  is  then  the  total  heat  output  of  the 
launch  vehicle  from  engine  Ignition  until  the  time  the  launch  vehicle 
passes  through  the  layer.  The  heat  output  has  already  been  adjusted  for 
the  effects  of  radiation,  after-burning,  and  the  deluge  water. 

In  the  acid  drop  deposition  model  the  vertical  velocity  used  to  carry 
the  drops  up  into  the  cloud  are  computed  assuming  that  there  exists  a 
parabolic  velocity  profile  with  a  mean  vertical  velocity  equal  to  the  cloud 
rise  velocity  computed  in  each  layer.  The  height  to  which  drops  are 
carried  is  computed  by  Integrating  the  drop  net  velocity  (local  cloud 
vertical  velocity  minus  the  settling  velocity  for  each  drop  size  category) 
in  time  for  as  long  as  the  net  velocity  is  upward.  This  integration  is 
performed  for  points  at  the  center  of  the  cloud  and  at  the  edge  and  defines 
the  locus  of  heights  reached  by  each  drop  size  category.  Drops  are  assumed 
to  fall  clear  of  the  cloud  when  their  net  velocity  becomes  negative.  Drops 
falling  out  in  each  layer  are  assigned  to  that  layer  as  a  source  for  the 
deposition  model.  Since  the  deposition  model  output  depends  upon  the 
assumed  drop  size  distribution  and  the  only  drop  size  distributions 
available  were  measured  within  the  cloud  at  heights  of  700  m,  the  initial 
(ground  level)  distribution  was  computed  by  using  the  cloud  rise  model 
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and  the  drop  source  algorithm  to  match  the  distribution  at  the  measured 
height  of  700  m. 

3.  REEDM  OUTPUT 

Examples  of  REEDM  output  are  shown  in  these  last  view-graphs  where 
isopleths  of  acid  drop  deposition  are  shown  on  a  maps  of  KSC  and  VAFB. 
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Figure  1,  Schematic  diagram  illustrating  major  components  of, the  REED 
structure.  “  —  . 
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FIGURE  A-1.  REEDM  exanple  run  meteorological  data  file,  taken  from  rawinsonde 
observation  at  Cape  Kennedy,  Florida  12  November  1981. 
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TABLE  2.  fuel  EXPENDITURE  AND  HEAT  CONTENT  DATA 


Vehicle  Type 

Property 

Space 

Shuttle 

Titan 

III 

Delta 

2914 

Delta 

3914 

(a)  Nonnal  Launch 

Fuel  Expenditure 
Rate  W  (g  s  ) 

1.5219x10^ 

5.437x10^ 

8.3607x10^ 

1.0576x10^ 

Effective  Fuel 
Heat  Content  H 
(cal  g  ) 

1479.7 

2021.1 

1766.0 

1449.9 

(b)  Launch  Failure 

Fuel  Expenditure 
Rate  W  (g  s  ) 

9.8873x10^ 

1.3594x10^ 

2.7294x10^ 

3.7073x10^ 

Effective  Fuel 
Heat  Content  H 
(cal  g  ) 

1000.0 

1000.0 

690.0 

411.2 

Burn  Time 
tfi  (s) 

1027.0 

240.0 

69.0 

126.0 

TABLE  3 .  VALUES  OF  THE  CONSTANTS  IN  THE  EXPRESSION  FOR 
VEHICLE  ALTITUDE  VERSUS  TIME 


Constant 


Vehicle  Type 


Space  Shuttle 
Titan  III 
Delta  2914 
Delta  3914 


0.652213 

0.429580 

0.922156 

1.245756 


0.468085 

0.518422 

0.432703 

0.418095 


0.375 

5.0 

0.54 

0 
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TABLE  4  .  EXHAUST  CLOUD  CONSTITUENTS  (FRACTION  BY  WEIGHT) 


Vehicle  Type 


Constituent 

Space 

Shuttle 

Titan 

III 

Delta 

2914 

Delta 

3914 

HCi 

0.0379 

0.1932 

0.1218 

0.1589 

A^2°3 

0.1828 

0.2819 

0.2214 

0.1936 

CO^ 

0.2503 

0.2665 

0.2055 

0.2783 

CO 


0.00042 


0.0222 


0.0156 


0.0331 
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stabilization  Time  for  STS-2 


Ground-level  HCA  gravitational  deposition  Isopleth  plot  using  the  "LAND"  map  of  Kennedy 
Space  Center, 
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